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The time evolution of the adiabatic piston problem and the consequences of its stochastic motion 
are investigated. The model is a one dimensional piston of mass M separating two ideal fluids 
made of point particles with mass m <^ M. For infinite systems it is shown that the piston evolves 
very rapidly toward a stationary nonequilibrium state with non zero average velocity even if the 
pressures are equal but the temperatures different on both sides of the piston. For finite system it 
is shown that the evolution takes place in two stages: first the system evolves rather rapidly and 
adiabatically toward a metastable state where the pressures are equal but the temperatures different; 
then the evolution proceeds extremely slowly toward the equilibrium state where both the pressures 
and the temperatures are equal. Numerical simulations of the model are presented. The results of 
the microscopical approach, the thermodynamical equations and the simulations are shown to be 
qualitatively in good agreement. 
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I. INTRODUCTION 

The adiabatic piston problem is a well-known example in thermodynamics which has given rise to continuous 
controversy for the last 40 years (see for selected references). It was also mentioned by E. Lieb as one of the 
"problems in statistical mechanics that I would like to see solved" at the StatPhys 20 meeting in 1998 [||. The 
problem is the following. The system is a finite cylinder containing two gases separated by an adiabatic movable 
piston. Initially, the piston is rigidly fixed by a brake and the two gases are in thermal equilibrium characterized by 
([/[I], V[i], A^[i]) and {U[2], V[2], N[2]) (see Fig. |l|). At a certain time to, the brake is released and the question is to find 
the final equilibrium state. 

Within the framework of thermostatics, one only knows that the final equilibrium state (oo), Vjij (oo), iV[i] (oo)) , 
(U[2]{oo),V[2]{oo), N[2]{oo)) , corresponds to a maximum of the entropy function S{U,V) = S[i] + S[2] = S[i]{U,V) + 
S[2]{Uo — U,V — Vb), where Uo = U[i] + U[2], Vq = V[i] + V[2], submitted to the constraints S[i] > S[i]{to), S[2] > 
S[2]{to). The solution of this mathematical problem yields the result that the final pressures must necessarily be equal 
(oo) = p[2] (oo) (j'-G- mechanical equilibrium), but nothing can be said about the final temperatures T[x] (oo), Tpj (oo) 
(originally Kubo had arrived at the wrong conclusion (oo)/T[i] (oo) — p[2] (oo)/T'[2] (oo) f^). In other words the laws 
of thermostatics are not sufficient to predict the final equilibrium state. Here started the controversy, because some 
physicists could not accept this limitation of thermostatics. Callen writes in |^ "the movable adiabatic wall presents 
a unique problem with subtleties" . The problem is however neither very subtle, nor unique since it appears every 
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FIG. 1. The adiabatic piston setting. 
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time movable adiabatic walls are present, and they are indeed essential in the axiomatic formulation of thermostatics 
(see e.g. g). 

Recently, with the help of an oversimplified macroscopical model for the gases, it was shown that the adiabatic 
piston problem, can be solved within the framework of thermodynamics 0. Equations for the time evolution were 
derived from the first and second law. They show that the equilibrium state depends in an essential manner on 
the viscosity of the gases and one has to solve explicitly the time evolution to find the final state, which is however 
uniquely defined. In particular in the final state the temperatures T[i](oo) and T'[2](oo) are different. 

Using arguments from the kinetic theory of gases, similar equations were previously obtained by B. Crosignani et 
al. 1^, without the use of the laws of thermodynamics, but neglecting the effects of fluctuations. The conclusion were 
however the same as the one obtained from thermodynamics, namely T[i](oo) 7^ T[2](cxd). At this point the controversy 
started again: it was argued that the stochastic motion of the piston, induced by the stochastic forces exerted by the 
particles of the gases, might lead to a "true" equilibrium state where T[i](oo) ~ r[2](oo) andp[i](oo) —p[2]{oo)- In ^ 
it was stated "this stochastic motion leads to a violation of the second law" and "this would be equivalent to permit 
heat to flow although the piston itself remains an insulator". It was also suggested that this evolution toward the 
"true" equilibrium state, if it exists, might be so slow that it can never be observed in practice. 

The work done by J. L. Lebowitz in 1959 should have given some insight into these questions already a long 
time ago. Indeed, using a model of non-interacting particles (in one dimension) to describe the gases, J. L. Lebowitz 
computed for the first time in |7j , the heat conductivity of a movable "adiabatic" piston at lowest order in e = -y/m/M 
where m is the mass of the gas particles and M is the mass of the piston. However, we should remark that the model 
discussed in [Q was slightly different than the one considered here since the piston was constrained to remain in 
mechanical equilibrium by means of an external force. 

This same microscopical model, but without external force, was reintroduced by Piasecki and Gruber in g. Starting 
with the Boltzmann equation, they obtained the exact stationary solution for the infinite one dimensional system, with 
= P[i] = P[2] = P[2]it)^ T[^{t) = T[i] ^ T[2] = T[2](i), in the very special case M = m. The main conclusion in ||] 
was that the space asymmetry {T]^i] ^ T[2] ) conspires with the stochastic motion to give a stationary non-equilibrium 
state with non-zero average velocity of the "piston" toward the region of the higher temperature. The same model 
with M = m was then considered by Piasecki and Sinai to study the finite system. They were able to show 
that the piston evolves toward an equilibrium position corresponding to uniform density and uniform distribution of 
velocities (thus uniform "temperature"), for the particles throughout the system. Simultaneously Gruber and Piasecki 
investigated the more physical situation m ^ M ]To| . Assuming that the solution of the Boltzmann equation has an 
asymptotic expansion in e = m/M , the stationary solution of the infinite system was computed at the order e^. 
It was again established that the temperature difference (T[x] 7^ ^[2]) induces a macroscopic motion toward the high 
temperature region despite the absence of any macroscopical force {jp[i\ = P[2])- In particular, it was shown that in 
this stationary state the piston behaves as a randomly moving particle with kinetic energy kB^JT\^i\T\^2]/'^^ average 

velocity ^y^(y7]^-y7]^). 

In all these recent articles, the physical problems of heat conduction, viscous force, work power, as well as the 
approach toward the stationary, or the equilibrium, state were not considered. These questions are investigated in 
the present article. 

Although we did not solve the problem mentioned by E. Lieb, which would need a rigorous mechanical many-body 
approach, the following analysis answers some of the conjectures discussed above and should contribute to a better 
understanding of the physical mechanisms involved in the time evolution of the adiabatic piston. It raises also some 
conjectures that one might be able to prove. 

The main conclusion of this paper is that the phenomenological macroscopic thermodynamics, the microscopical 
Boltzmann equation, and the numerical simulations, give results which coincide with a very high degree of accuracy 
when m <C M . In other words there is no violation of the second law and no paradoxes involved. 

The results of this investigation show that a wall which is "adiabatic" when rigidly fixed (no heat transfer), 
becomes "conducting" as soon as it is allowed to have a stochastic motion, and this is independent of the fact 
that its macroscopic velocity is zero or non-zero. This means that in this model heat transfer is associated with the 
stochastic motion of the piston, while work power is related to the macroscopic (or average) velocity of the piston 
{V)^. The time evolution obtained from thermodynamics (where the above conductivity is taken into account), from 
Boltzmann equation, and from the numerical simulations, shows that the system evolves toward a final state where 
the temperatures, as well as the pressures, are equal. However, and this is the fundamental point, when m <C M the 
time scale needed to reach this final state is so enormous (about the age of the universe for realistic numbers), that no 
reasonable person would call such a piston a conductor. More precisely our numerical analysis indicates that as soon as 
m/M is small enough, the time evolution proceeds in two stages with time scales which are very different. In the first 
stage the evolution is adiabatic and proceeds rather rapidly, with or without oscillations, towards a state of mechanical 
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equilibrium where the pressures become equal, but the temperatures remain different (this is the "quasi-equilibrium 
state" suggested in Q). After this quasi-equilibrium situation is reached the piston will slowly drift into the high 
temperature domain and the system evolves, on a time scale several orders of magnitude larger, with equal pressures 
and slowly varying temperatures toward the final equilibrium state where temperatures, densities and pressures are 
equal on both sides of the piston. On the other hand for m/M ~ 1, the adiabatic piston is a good conductor and 
thermal equilibrium T[i](oo) = T[2](oo) is reached rapidly. Notice that on a very short time scale, there exists a stage 
where the piston behaves as if the system is infinite. In this stage, if the initial pressures are different, the piston is 
accelerated while if the pressures are initially equal the piston evolves with a (small) constant velocity inside the high 
temperature fluid. 

In the following analysis the microscopical model is a one dimensional system of point particles which do not interact 
except for purely elastic collisions. It is exactly this model which has been simulated on a computer. 

For the analytical analysis we introduce three main assumptions. First we assume that the correlations between 
the velocities of the fluid particles and the velocity of the piston can be neglected. Then we assume that the velocity 
distribution functions for the fluid particles on both faces of the piston admit a scaling function with respect to their 
mass m, as it is the case with Maxwellian distributions. However, because of the simulations (Sec. VIII ), and the exact 
solution of Piasecki and Sinai we do not assume Maxwellian distribution for the velocity of the fluid particles. 
The temperatures T^ij {t) and Tp] (0 '^'^ both faces of the piston are then defined by the second moment of the velocity 
of the fluid particles. Finally we assume that the moments (y^)^ of the velocity V of the piston, s > 1, have an 
asymptotic expansion in powers of e = yJm/M . 

To discuss the solutions of the time evolution equations, we shall either assume that the fluids densities and 
temperatures are constant (in the case of an infinite syst em) , or uniform (in the case of a finite system). 

The microscopical model is defined in Sec.||. In Sec. Ill we use the conservation laws of the problem to derive 
the dynamical equations for the first and second moment of the piston velocity distribution and give in Sec. IV 
these equations for arbitrary moments (V^'')^. The relation between the microscopical model and a phenomenological 
thermodynamic approach is discussed in Sec.^ where we derive a microscopical definition of the friction coefficients 
and the heat conductivity. The analysis of the stationary non-equilibrium, and the equilibrium, states is presented in 
Sec. VI, while the time evolution is discussed in Sec. VII, The numer ical s imulations, the comparison with the d yna mic 
equations, and the thermodynamic equations, is the subject of Sec. VIII. Finally, conclusions are given in Sec. IX. 



II. MICROSCOPICAL MODEL 



We consider the fluids to be perfect gases made of A^[i] , respectively N\^2\ , identical non-interacting point particles 
of mass m, making purely elastic collisions, and contained in a flxed cylinder of length L. The basis of the cylinder is 
a circle with area A. The adiabatic piston is a rigid solid (i.e. no internal degree of freedom), with mass M, described 
by the microscopical variables [X, V = ^) (see Fig. |l|). 

We assume that the cylinder is an ideal constraint, which has no physical properties: the collisions of the particles 
on the cylinder are purely elastic, and there is no friction between the piston and the cylinder. 

For t < 0, the piston is flxed {X = Xq, V = Q) and both fluids are in thermal equilibrium, with temperatures Tjij, 
respectively T^j, described by the Maxwellian distribution functions: 

P[i](a:,w;i < 0) = p[i](/)T[i, (w) for < x < 

P[2](a:,w;i < 0) P[2](/>Tp,(w) for Xa<x<L (1) 



where 



and 



l^[l] 1 1JV[2] 1 



At time t = 0, the break is released and the piston can move freely. It will thus have a stochastic motion induced 
by collisions with the fluids and (X, V) become random variables. The problem is to investigate the random motion 
of the piston. To simplify the discussion, we consider the system to be one-dimensional. In this case the adiabatic 
piston is just a point particle with mass M . 
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For t > 0, the piston moves freely. Since the collisions are assumed to be purely elastic, then if the masses M and 
m have velocities V and v before the collision, their velocities after the collision are given by 

2M ^ 
V =v — — [v — V ) 

1m 

V' = V + — {v-V) (4) 

The main assumption of our analysis is that the time evolution is such that it is possible to neglect the correlations 
between the velocity of the piston and the velocities of the fluids particles. This means that we assume that the joint 
probability distribution to find, at time the piston at position X with velocity V and a fluid particle at position 
X ± |e| with velocity v can be expressed by 

<f>[i] (X, F, X, t) = t/; t) Imi {X - |e| , v; t) = <i>(X, V- t)n[i] (v, t) 

a>[2] [X, V, X, V- 1) = V- 1) Imi p[2] {X + |e| , v; t) = ^{X, V- t)n[2] (i)0[2] {v. t) (5) 

where we have introduced the time dependent particle densities on the faces of the piston njjj and the normalized 
fluids velocity distribution functions on the faces of the piston 0[j](t),t) for i = 1,2. 

In the case of an infinite cylinder, we shall consider that and «[2](i) are constant, equal to the fluid initial 

densities, and that = (pT^^iv) {i = 1,2). This assumption is justified if m <C M, since in this case there is 

a vanishing probability that the piston interact back with the perturbation it causes in the state of the surrounding 
fluids (no recoUision). In other words, the piston always "sees" on both sides the unperturbed initial states. 

For the finite system, when m <C M, we expect the average velocity of the piston to be small, and therefore the 
time needed for a macroscopical but small displacement of the piston will be much larger than the relaxation time of 



the fluids (if the initial pressures difference is not too large). We shall then assume in Sec. IV that at all time t, the 
fluids on both side of the piston are homogeneous and thus that the distribution functions p[^(x,v,t), < x < X, 
and p<^2]{x,v,t), X < x < L, are independent from the position x. We then have 

«w(t) = ^, = (6) 

and the time-dependent temperatures Tji] (t) and T[2] (t) a-re defined by the second moment of the fluid velocity. 
For the following analysis, we introduce the small parameters 

2m , [rn 

« = TT and e — \—. (7) 

M + m y M ^ ^ 



III. CONSERVATION LAWS 



We use the conservation laws present in the problem to determine the dynamical equations for moments of the 
piston velocity distribution. In particular, the conservation of momentum and energy leads to equations for the first 
and second moment of the piston velocity distribution. 



A. Linear momentum 



Let (t) and n[2] (t) denote the linear momentum of the fluids on the left and on the right of the piston at time 
t. Using the collision equations (||) and the assumption (^) the force exerted by the left fluid on the piston is 

F^^^^P = - \im(—^) ^ Ma J dV^,{V;t) dv n[,]{t)4,[^{v,t)(v - V)^ (8) 

where ^^{V;t) is the time dependent velocity distribution function of the piston. Similarly the force exerted by the 
right fluid on the piston is 

F[2Hp = - hm / J^\ = -May dV<i>,{V;t) J dv n[2]m[2]{v,t){v - V)^ (9) 
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Therefore, the conservation law of hnear momentum 



dt 



can be written in the foUowing form 



M— (V)f + hm 



6t 



St 



Using the expansion 



-{V)^^al dV<^>,{V;t) I d« [n[i](i)0[i](t;,t)(t;-y)2-n[2](t)0[2](-t;,t)(i; + y)2 



a / dV<^>e{V; t) / dv{v~ Vf [n^^ [t)c^[^ {v, t) + np] {t)c^[2] {v, t)] 



where 



(^ = 1,2), 



v=0 



we obtain 



-ly),=a[ny,^{t){v')\'^-ny,^{t){vy^ 

~{V),2a [7i[i](t) («)W+n[2](i) (t;) 

+ (y^)^a[n,](t)(.0)W-n„W(.°)f 



2a 



fc=0 



(fc + 3)! 



n[i]W0;J)(i)+np](t)<^;^]^(t) 



'(*=)/ 



where for any non negative integers s 



dv<j)[i]{v,t)v'', 
dv 4>\2\{-v,t)v'' 

/oo 
dV<^e{V;t)V' 
-OC 



and 



(10) 



(11) 



(12) 



(13) 



(14) 



(15) 



We should insist on the fact that in the definitions of the integration is on v positive. In particular (f'^)''' is 

not equal to 1. 

We should note that the evolution equation (^^ for the average velocity of the piston depends on higher moments 
of the piston velocity. We should also remark that if the term n[i](i) — 'n-[2]{t) is zero (which means, 

as we shall see in Sec.^ that the pressures exerted by the fluids are equal) and if we neglect the fluctuations of the 
piston velocity ((V')^ = {V)^) we find the stable equilibrium solution {V)^ = 0. The effects of the fluctuations is the 
main problem to be investigated and will give a very different picture. 



B. Energy 

Similarly, let (t) and E[2] {t) denote the energy of the fluids on the left and on the right of the piston at time 
t. It is equal to the kinetic energy, since the particles do not interact except for elastic collisions. Using the collision 
equations ^ and the assumption (||), we have 
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<5t-.0 \ St 



and 



lim 



6t 



dVMV;t) 1^ dvn[^{t)(l)[^{v,t){v ~Vy [V + -{v ~V) 



V + -{v-V) 



Therefore, the conservation of energy 



Ma I dV^,{V-t) I dvn^2]m[2]{v.t){v~Vf 

— OQ 



2 dt^ st^a 



St 



St 



= 



leads to 



dt 



(V^)^ = 2a dV<i>e{V;t) / dv 



{n[i]im[i]{v,t)iv-Vf 



V ■ 



/ dv{v^Vf[V+-{v^V)\ hi](t)0[i](,;,i)+n[2](i)'/'[2](^^,i)] . 



-2a 

With the expansion ( |12| ) and the definitions (|f5|) this last expression can be written as 
d 



dt ^ 



n,,,{i^{v^)f+n,2,{t){v^)\ 



+ (y), a{2 - 3a) [n^] (i) (t;^); ' - n^, (i) (t,^) 
- (1-2)^ a(4 - 3a) (t) (t;) « + np] (i) (i^) 



[2] 



+ (F=^)^a(2-a) n^,^{t){v^Yl^-n^2]{t){v') 



2a 



fc=0 



(fc + 4)! 



(2fc + 8 - 3a) n[i] (O^jir (0 + np] (t)4, ' (t) 



(16) 



(17) 



(18) 



(19) 



(20) 



We thus obtained an equation for the evolution of the second moment of the piston velocity. Let us remark that 
even in case of mechanical equilibrium where n^i^it) — "-[2](i) = the initial state (V^*)o = is unstable 

if ?T.[i](i) (v'^)^^^ + iT'[2]{t) 7^ 0- conclusion, if we do not neglect the correlations in (|l^), we see that even in 

case of mechanical equilibrium, the initial state (^^^)o = is unstable and thus the piston will start to move. 



IV. EQUATION FOR ARBITRARY MOMENT 

In the previous section, we derived equations for the first and second moments of the piston velocity which are 
coupled with higher order moments. We derive in this section the infinite set of coupled equations for all the moments 
of the piston velocity. 

We use the Boltzmann's equation for the piston velocity distribution derived in ||] under the assumption (|^) 
-g^'^eiV; t)^ j dv\V^ v\ [Q{V - i;)n[i] (t)0[i] {v-{2- a){v - V)- 1) {V + a{v -V)-t) 



+Q{v ~ V)n^2] (i)0[2] - (2 - a){v -V);t) $e {V + a{v ~V);t) 
^eiv - T/)n[i] (t)0[i] («, t)$,(V^, t)-e{V- «)np] (i)'^[2] {v, ^eiV, t)] 



(21) 



to obtain the equations for the velocity moments 
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dvn[^{t)(t>[^{v,t)iv - V) {[V{1 -a)+ av]' - V'} 



V 

-l^j dvn^2]{m[2]{vM^-V){[V{l-a)+avY -V'}^ (22) 

where {A{V))^ = ^^^dV A{V)<^>e{V-t). This equation can be written in the foUowing form 

I {V')t - dvn^,^{t)cj>[^^{v,t){v - V) {[V{1 - «) + avr ~ V^}^ 

+ l^j^ dv n[2] (i)0[2] {-V, tKv + V) {[V{1 - a) - av]' - V'}'^ 

dv {v ^V){[V{l-a)+ avY - V'} [n[i] (t)</)[i] {v, t) + n^^] ^[2] {v, t)] ) . (23) 



Using the expansions (|lj) and 



[V + aiv - V)r - - E - VyV^-", (24) 



we obtain 



-{V^)^ = Y^AsAt){Vn, (25) 

r=0 

where one has the expUcit expression of As,rit) as a function of the parameter a. In the following, it will be more 
convenient to have the explicit expression of As^r{t) as a function of the parameter e = y^m/M. 

Motivated by the Maxwellian distribution (||), wc introduce the scaled velocity w = ^frnv and the corresponding 
normalized distribution functions f[i](w,t) 



(w, t) = ^/m/[,](^/mz;,t), / dw f^i^{w,t) ^ 1, (i = l,2). (26) 

'J —OO 

We then have from Eqs. ( p^ , [2^ ), 

poo 

W = m-«/2 / dwf[^ {w; t)w'i = m^'/^ (^9) W ^ ^ 2) (27) 
Ja 

and 

^(^(O^me^+i^/Vff (t). (28) 

For any integer r we introduce the functions Kr\t) and Kr(t) defined in the following manner. For any q > 0, and 
i = l,2, 

K%,{t) = M-^/\,^it){w^)^\ (29) 

if[!U^(t)^M(^+i)/2n[,](t)/[y)(t), (30) 

if,(<) = xW(<)-(-l)''ifPl(t) forr>-2, (31) 
Kr (t) = X^^l (<) + ifl^l (t) for r < -3. (32) 

With these definitions, we have 



7 



, cso 

-(n* = E^«.^w(^'')* (33) 



dt 

r=0 

where the coefficients Ag.rit) are given by 

AsAt) = ^^^^PsA^^)Ks-r-i{t). (34) 

The functions Kr{t), Eqs. ( |3l| , p2| ), depend only on the densities and the velocity distribution functions of the fluid 
particles on both faces of the piston, while Pg^ri^^) is a polynomial in e^, of order min(r, s — 1), with constant 
coefficients, given by 

P,Ae') - (1 - e'y-' [{s + 1 - r) - e\s + 1 + r)] ,. ^^""''^ r < . - 1; (35) 

■' r!(s — r + 1)! 

(37) 

P......(.^) = -2g(-l)n2.Y(l + -^)-- \._i:'^y,'],^3)! ^ ^ + 



One can verify that Eq.(|33|) for s = 1,2 are identical to the equations for the first and the second moment (|14|,|20|) 
derived in the preceding section. 

V. FROM KINETIC THEORY TO THERMODYNAMICS: FRICTION COEFFICIENT AND HEAT 

CONDUCTION 

In this section we relate the microscopical model to a phenomenological thermodynamic approach. This will suggest 
microscopical definitions for "pressure" , "viscous force" , "temperature" and "heat flux" . 

Using a very primitive model for the fluids, the equations for the time evolution of the piston problem were derived 
in 0] from the first and second laws of the thermodynamics. They read 

(A/ + <5A/)^=A(p[i]-p[2])-(A[i]+A[2])T/; V=^ (39) 



dS[i] _ A[i] 



K 



V' + 7^(T[2]-r[i]) (40) 



dt T[i] T[i 

dS[2] A[2] 2 K 

where X and V are the macroscopic position and velocity of the piston of mass M, SMV'^/2 represents the kinetic 
energy of the fluids, A is the area of the cylinder, P[i],P[2], A[i], A[2], k are phenomenological time-dependent functions 
of the variables (X, V, »S'[2]) and denote the pressures, the friction coefhcients, and the heat conductivity. S[i] and 
S[2] denote the entropy of the fluids on the left and on the right of the piston. 

For ideal fluids in one dimension, the thermal equation ( pO[ ) and (^) are equivalent to 

l^inkB^^ - hnV' - + .(T,2, - T,,) (42) 

. A„ + ^M^y „ .(T,, _ T,,). (43) 

In thermodynamics, the adiabatic piston is defined by k = 0. One should remark that in this case (k = 0), if 
the velocity of the piston is sufficiently small so that the term in can be neglected, the evolution of the piston is 
damped, but reversible in the thermodynamical sense, characterized by the adiabats 
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= 5'[i] and S'[2](<) = 



j.e. 



T[^X^ = etc 



T[2](L-X)2 = cte. 



(44) 



In particular, the solution = p[2], = is stable. 

To make the connection with the thermodynamical equations (|39|), we write the microscopic equation for conser- 
vation of momentum ( |ll|) or ( |l^ in the form 



-(F),2mn[i](i) 



1 



(F),2mn[2](t) 



2(.)fi + (y),(«")f 



1 



(1^>* Jo 



fc+3 
t 



4™ 



fe=0 



(fc + 3)! 



"[i](i)0S'(i)+^[2](i)0;J^(i) 



(45) 



We are thus led to identify {V)^ with the macroscopic velocity V{t) of the piston, and to define the "pressure", the 
"temperature", and the "friction coefficient", on both sides {i = 1,2) of the piston by 



P[,](i) = n[,](t)fcsr[,](t) = 2mn[,](i) {v^)f = 2n^^{t) {wYt = 2MK^\t) 



(46) 



= 2mn[i](i) 

which gives in particular 



A[,](T/ = 0) = 4V^nH = 4eAfK[!\, 



(47) 



(48) 



A[i](V") + X[2]{V) = M [AeK^i - 2e^K^2V + 0(6^^2))] 



(49) 



and shows that the condition A[i](l^) > will be satisfied in the domain considered here where the average piston 

velocity is much smaller than the average absolute velocity of a fluid particle, V <C {v)^^^. Notice that the temperature 
definition Eq.(^) is consistent with a Maxwellian definition of the temperature. We thus have from Eq. (p5|,p|,p 10) 



M{l + e')-{V), = (p[i](i)-p[2](<)) - (A[i](y,t) + A[2](y,t)) {V), + F,,{t), 



and 



FW^P(i) = -FP^W(i) = 



1 



l + e2 
1 

l + e2 



[1] 



P[i](t)-A[i](T/,0 {V), + Fl 
-P[2]{t)-X[2]{V,t){V),+F, 



[2] 



where the stochastic forces are 



[1] _ 



fc=0 



[2] 



oo 

{V^\ (V)^] 2mn„(t) - ^ [(^'^+^). " 



/c+3 
t 



(fc + 3)! 
4m 



fc=0 



(fc + 3) 



■"[2](i)0;j^(o 



and^^st^Fit'+^if 



(50) 

(51) 
(52) 



(53) 
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We now compute the heat flux from one side of the piston to the other side in the microscopic model. This 
wiU give us a microscopic definition of the heat conductivity n. We wih see that the microscopical equations are 
compatible with the phenomenological thermodynamical equations in which we take into account a heat flux charac- 
terized by a conductivity k 7^ which comes from the fluctuations of the microscopical model. From the flrst law of 
thermodynamics 



pP^[l] , p[2]^[l] 



(54) 



where P^^'^\ respectively P^'^'^', denote the power transmitted to the fluid on the left in the form of work, resp. 
in the form of heat, together with the thermodynamic deflnition of the work power, i.e. 

P^^^\t)^FP-^^\t)V{t) (55) 

we are led to deflne, through Eqs. (|5l| , |l6| , |l8| , ^3|) , the microscopic heat power, or heat flux, and the microscopic heat 
conductivity k 



P, 



[2] -►[I] 



Q 



{t) = 



dE[i] 
dt 



1 °° r 1 



K 



[1] 



Tr 



(56) 



where aI^I are defined by Eq.(||) with K^l^. 



At order 0(e), it reduces to 



(57) 



To compare our definition ( |5^ ) with the heat conductivity obtained in 0| , let us consider the case of the infinite 

system. In this case the coefficients A'r ' and the densities n[i] are independent on time. As we shall see in Sec.0 for 
any stationary state of the inflnite system we have 



2K-1 



0{e) 



which thus gives for a stationary state 



□[2]^[i] ks 8mkB 



In particular, for a Maxwellian distribution of the fluid particles, this last equation reads 

= (T[2]-r[i])+0(e2). 

'nil v-' 111 ^ "12] 12] 

Therefore the heat conductivity is 

.„ Is^u^ V^[i]^[2] i-0(e2) 



M 



K = —\l "[1]"[2] 



"WV^I +"[2]\A[2] 



(58) 



(59) 



(60) 



(61) 
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Similarly, for the infinite system the friction coefficient (Evh can be written as 



\i]iV) = n 



[1] 



TniTO 



/ 2m 1 



dv exp 



mv 



(62) 



It is interesting to note that for a stationary state of the infinite system and Maxwellian distribution for the fluids, 
the expression for the friction coefficient X{V = 0) Eq. (|6^), and for the heat conductivity k Eq. (^Tj) are, to lowest 
order in e, identical with those first derived by J..L. Lebowitz for a slightly different model |^]. In 0| the piston 
is constrained in a state of mechanical equilibrium (i.e. {V'^} — for all odd integer k) by mean of an external 
potential and the heat flux is just given by dE^^/dt. In our model there is no external potential and in the stationary 
state the piston has a constant non-zero average velocity; we thus have to take into account the work power in the 
definition of the heat flux Eq. (|5^ . We should also remark that in j7| the friction coefficient is introduced via the 
Ornstein-Uhlenbeck process describing the motion of the piston, while our definition is associated with the mechanical 
equation of motion Eq.(^9|) and thus the factor M^^ difference between A in and our \{V — 0), Eq.(|6|). Finally, 
we want to stress that in the definition ( ^7| ) and ( |5^ ) of A and k, the fluid particles can have any arbitrary distribution 
and the system is not necessarily in a stationary state. 

To conclude this section we remark that we can use the thermodynamic equation 



dE\ 



[1] 



dX dXd^X 

dt + 

to define the time evolution of the microscopic entropy by 

dSti 



(63) 



Tr 



[1]- 



dt 



I dE,,] \ d 



M 



-- E {Vlt + {V), - m {V), - {V) 



k=0 



M 



^ll] 



[2(«;)WM(0^-(^3^W]+0(a 



(64) 



i.e. with (13) and (E^) 



dt 



m 
M 



2(«;)W 



{V)t+Pc 



Q ' 



(65) 



which is consistent with the thermodynamical equation (pO[). Again, assuming a Maxwellian distribution for the fluid 
particles, we will have 



Tr 



In summary, if initially (V^'^)j^q — for all r > 1, then as long as the contributions 



M 



M{{V^)^^{V)l)-kBT[^ , (z = l,2). 



(66) 



(67) 



as well as the stochastic force Fgt, can be neglected, the time evolution Eqs.(^,^2) ^^^^ coincide with the thermo- 
dynamic evolution of the adiabatic piston Eq.( p9|^d|) with k = 0. This will be the case as long as {V^)^ ~ 
and » fcBmax(r[i],T[2])/M. Furthermore, as long as {V)f ^ 1 the evolution is described by the adiabats 

S[i] (t) — S[i] , S[2] (t) = S[2] ■ Therefore during this first stage the piston evolves according to ( |3^ ) towards a state 
of mechanical equilibrium, i.e. p[i] ~ p[2], but not thermal equilibrium [T^i] ^ '^[2])- As soon as the pressures are 
approximately equal, the above assumptions are no longer valid, and as we shall see in Sec. VI] the system will evolve 
toward a state of thermal equilibrium. Numerical simulations will be presented in Sec. VIII to show the time scales 
on which these evolutions take place. 
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VI. STATIONARY NON-EQUILIBRIUM AND EQUILIBRIUM STATE 



A. Preliminary remarks 



We have seen in Sec. [V that the time evohition of the piston is characterized by the equations 

, oo 
r— 1 



where Cs{t) = ^s,o(^) a-nd for any r > 



e 



\s-r-l\ 



In particular for s = 1 and s ~ 2 we have (where = denotes "for MaxweUian distributions"), 



with 



P[i] ~P[2] = M2Ka = fcs(n[i]T[i] - n[2]T[2]) 



0)\# fc- * ^ /8fcs, / ^3/2 , „3/2 

2K^2 = "[1] - ?^[2] 



To lowest order, i.e. in the limit e — > 0, Eq.(68) together with the definition ( [46| ) of the pressure yields 

As discussed in Sec.||, for an infinite system with e ^ 1, the assumptions 

n-li] (t) = n^i] , (i;, t) = (/)T[., (v), (i = 1, 2) 

are justified and give (t) = , p[2] (t) — pp] constant. 
In this case, for the initial conditions 

we obtain at lowest order 

and thus there is no stationary state in the infinite system unless — p[2] in the limit e ^ 0. 

On the other hand, for the finite system, if the evolution is sufficiently slow and e ^ 1, we may assume 
following homogeneity conditions hold: 

""W^^^^W;' "W(^) = ^(1)^ ' (i?w)^ = %m(i;2)W EEiiV[,]fcBr[,](i), (z = l,2) 
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where we do not assume Maxwellian distribution of the velocities but define the temperatures by the second moments 
of the fluid velocity distribution Eq.(|4^). Under this homogeneity condition, the time evolution for T[i] and T[2] are 
therefore 



dt 



N[i\kB 



M 



(78) 



In the limit e 



^^^^)t =kB 



we have from Eqs.(||,0,^ 

dt 

dT,2] 



3}i 



Tr 



[1] 



L 

d 



-Tr 



[2] 



dt 



L-{X) 



which yields the "adiabats" 



T[i](t) {Xf, - 7[i] = cte; r[2](t)(L - {X),f = 7[2] 
and the undamped oscillatory motion 



cte 



2]7[2] 



iL-{X),r 



(79) 



(80) 



(81) 



In conclusion, to lowest order in e, there is no approach toward an equilibrium state in the finite system while a 
stationary state in the infinite system can exist only if the initial pressures are equal (p[i] = P[2])- This conclusion is 
consistent with the previous result, since for finite M the condition e — > means that to ^ and thus the friction 
coefficient and the heat conductivity are both zero. 



B. Stationary state of the infinite system 



In this section, we discuss the stationary states of the infinite system. In this case the densities n[2], and the 
coefficients Kr are constant. Moreover we consider the moments {w'^) which appear in the definition of Kr to be given 
by the Maxwellian distributions i.e. 



1 



1 



(-°) =i> (-^) =^v2^bt,], 

= ifc^T[,], (^.3)H^^(2fc^J.^,^])3/^ (^ = l,2). 
The stationary states of the infinite system are given by the solution of the equations 

o = ^ei-'-iiF,,.(6^)i^,_._i(noo- 

From Eqs.( [70| , [7ll) , it follows immediately that 



X 



1 - (1 - 2e ) ' \ --(^[1] - n[2]. 



(1 - 2e2 



P\l]-P[2] 



2-MKi 
A 



i.e. at lowest order in e we have 

P[l]-P[2] - P[2] \ ^ TO 1 (P[l] P[2]\ m P[l\Tl2] - P[2]T\i] 



A 

- Pm 



A 



M 



A fcs Vr[i] r[2]/ M A^7]I^ 



(82) 
(83) 

(84) 
(85) 

(86) 
(87) 
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which shows that is impossible to have a stationary state with (V^);^ = 0, = p[2], T^i] ^ Tp]. 
In the following we shall assume that the solution of Eq. ( ^3|) has an asymptotic expansion in e 



;=o 



Then in the limit e — s- we must have 



which can also be seen explicitly with Eqs. 



Ko = 0, i.e. py 



\(0 

/ oo 



'■P\2] 



(89) 



since X = 0{e). 

To take into account the possibility of non-zero pressure difference for the infinite system, and to have results valid 
for finite systems where the equilibrium densities might depend on e, we now consider that the densities n[2] are 
function of e 



which then implies that the coefhcients Kr are also functions of e 

oo 



/=0 



The condition ( |8^ ) for a stationary solution is now 

'(0) 



K'o' = lim/!:o(e) =0 



I.e. 



"[1] ^[1] - "[2] ^[2] = 



= 



lim[p[i](e) -p[2](e)] = pf^ - p^^^ = 0. 
Therefore Eq.(^) can be simplified by e to give 



= p. 



s.s—1 
s-2 



oo 



s,2s — s— r 



We then solve ( |95| ) at successive order in e. 
At order e*^, we have 



= 4if(i) (F)^::) + 4K["^ 8Ki") (F^)^) (. = 2) 

= 2s4'^ (F-i)^^ + 2.(. - ) (F-^)^^ - 4.xi°) (. > 3) 



which yields 



'[1] V'"/ ^ "[2] 

P[i](e) -P[2](e) 



—J [11 J [91 hm 



" e (e) a/7[2] + P[2] (e) 



(90) 

(91) 

(92) 

(93) 
(94) 



(95) 



(96) 



(97) 
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K 



(0) 



1 



M 



(0) 



[2] 



T[i]T[2] lim 



^0 P[i] (e) + Pm (^) 



and 



=1 



[s/2] 

\(0) = V . 

loo 2''k\(s-2k)\ ° 



Therefore at order e'^, the distribution function for the velocity of the piston is 



1 



\/27rA 



exp 



2A 



Let us remark that we recover Eqs.( p6| , p7| ) at order zero. 
At order e^, introducing the notation 



1,(0 



(0 



„(i) 



2i^iT 



(0) 



Eq.(p5h yields 



-2fc«Vo + fcf +fci"^(A + yo') 



+(. - i)fc« (y-^);;' + |(,s - i)(. - 2)4") {v^-^y^ + k^^l 



s-2\(0) 
' oo 

As one can verify the distribution function at order is 

1 



(0) 

OO 



exp 



Tikf 



X {1 + e [a,{V - Vo) + a^iV^ - A - V„^) + a^iV^ - 3AVo - V,f)] } 
where Vq, A are given by Eqs.@^ and 



aiA3 = k^Yl2VoA^ - k['>VoA ~ ^fcf (A - V^^) 



SagA^ 



fci")A2 



„(o) 



(99) 



(100) 



(101) 

(102) 
(103) 

(104) 



Let us remark that for 



P[i] 



(105) 

(106) 
(107) 
(108) 
One can then 



P[2] then Vo = and we recover the results of Gruber and Piasecki 
successively obtain higher order terms for the velocity moments. 

In conclusion, for the infinite system we observe a stationary state with non-zero average velocity (V)^, which depends 
on n[2], T^i], T[2], given at order e by Eqs.(p6,102). 

It is interesting to note that this result present some analogy and differences with those obtained by J.L. Lebowitz 
in pi, where he considered an adiabatic piston constrained around the origin by an external force: The coefficient A 
Eq.(38|) is exactly the same as in 0; however for p[i] — p[2] = 0{e) ^ we have a drift Vb, while in Q the external 

force give {V^'"^'^)^^ = {k > 0). 

To conclude the discussion on the stationary states of the infinite system, we investigate under what condition there 
exists a stationary state with zero average velocity i.e. {V)^ = 0. 
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From Eq.(Q) it follows immediately that 

P[i] - P[2] = -ni{n[i] - "[2]) (V^)^^ + 0{e^). 

However one ha s to check that a stationary solution of Eq.(|8^) does exist with (V)^ = at all order in e. 
From Eqs.(||,^, {V)^^0 up to order e iff 



(0) _ 7^(1) _ 



K^^'^O, i.e. lim(p[i](e)-p[2](e))A = 



(109) 



(110) 



and 



(2) 



Expressed differently, (V)^ = up to order e iff 



(Ill) 



,(0)r 



"[1] ^[1] 



(l)r 



21 ^ [2] 



,(2) 



'[21 ^[2 



M 



(0) (0) 
"[1] "[2] 



(112) 



This last equation is equivalent to 



and, in this case 



P[l] - P[2] 



(0) _ (0) 
"[1] "-[2] 



(0) 

oo 



2\(1) _ 1.(1) 

oo 



At order e^, we have (taking s = 1) 



Therefore (y)^ = up to order iff 



K, 



(3) 



/f^A + ifi;2^(fci^^ - 2k'_'>A) = 



'(0)^1,(1) 



(1), 



Finally at order e'^, we have (taking s = 1) 



(113) 

(114) 
(115) 

(116) 

(117) 

(118) 



r>3 



Therefore (V)^ = up to order iff all K^^aie uniquely determined by the functions obtained at previous order. 
In conclusion, under the assumption Eq.(88) and given the density n[2], and the distributions (/)[i](w), (f>[2]{v), with 
7^ r[2], then the average velocity of the piston in the stationary state is non-zero for all densities n^^ except for 
one special value = "•[!](£) for which (V)^ = 0. For this special solution, one has 



(0) / 2\ [1] / 2\ 



2\[2] 



J_„('?) /,„2\[11 



(9>3). 



(119) 
(120) 
(121) 
(122) 



r>3 
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In this state we have (see Sec.^) 



1 , ,,[2] (w^fUwf^-(w^f^{wf^ 

and thus there is a constant heat flux for any distribution <j)[i]{v), (t>[2]{v) such that ^w'^)'"'^' (w)'^' ^ (^w"^)^^' {'w)^^\ For 
Maxwelhan distribution, this means that there is a constant flux for any T[i] ^ T[2] . 

Therefore there exists a special non-equihbrium stationary state with zero average velocity of the piston, which is 
however distinct from the one given in since in our case {V'^'^^'^)^ is not zero (for k > 1). 



C. Equilibrium state of the finite system 



Any equilibrium state of the flnite system is necessarily a stationary state with (V)^ ~ and p^l~*[^l = where, 
for (^)^ = (see Eq.(p|)), 



p[2]^[l] _ _ 



2M 



(l + e2)2 



Taking into account (|98|), Eq.(124) implies at order e 



/ 3\[ll / \[2] / 3\[2] / \ 



[1] 



(124) 



(125) 



Assuming Maxwellian distribution for the fluid particles, it implies r[i](oo) = r[2](oo) at order e and in this case we 
have 



(126) 



We conjecture that for any equilibrium state th e rela tion (126) must be satisfied, which then implies T^i] (oo) = T[2] (oo). 
Assuming this conjecture (126) to hold, Eq.(llO) yields first 



K^J\oo)^Q for 1^0,1 
i.e. njjj(oo) = nj2](oo) for ^ = 0, 1 



Then, with Eq.(|ll2|), Eq.(|l28|) is also vafid for / = 3, thus 

(oo) = for all r even, / = 0, 1, 2 

and with (|l0|) 

/y2fe+i\(i) ^ for all A: > 0. 

\ / OO — 



(127) 
(128) 

(129) 

(130) 



Iterating the argument at successive order we arrive at the following conclusion. Under the assumption ( pq ) and 
the conjecture ( |126D , any equilibrium state of the finite system must satisfy n[i](cx)) — n[2](oo), in other words 
P[i](oo) =P[2](oo)and T[i](oo) = T[2](oo). 

Conversely, if (v) — (f)^2] thus T^^ — T[2] , and n^^ (oo) = n^2] (oo) then the solution of the stationary equation 
at order e is 



(y2fe+i\, 







(y2fe^^ = (2A:- 1)!!A^ A = 



gi(oo) 
2ii:_i(oo) 



(131) 



Therefore at order e the velocity distribution of the piston is necessarily Maxwellian with temperature T — A/A/fc^, 
i.e. 
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27rA 



exp 



2A 



(132) 



In conclusion the adiabatic piston is stricto senso a conductor since the only equilibrium states are those for which 
T[i](oo) — T[2]{oo). However the important question to be considered in the next section concerns the time scale 
necessary to reach this true equilibrium state. 

VII. TIME EVOLUTION 

In this section we give a qualitative (non rigorous) discussion of the time evolution for the infinite and the finite 
systems. 



A. Infinite systems 



As discussed in Sec.||, for infinite systems and e <C 1 the assumptions 7t[,;](i) = ^^[i], (f>ii]{v, t) = (j)T^i] (v), i — 1,2 are 
justified. In this case the functions Kr are constant, and we are led to solve the equation (pSh 



where 



i.e. with t' = et 



-V^e[C + AV], V(< = 0) = 0, 



i[i] 



g|s-r-l| ^ 



(133) 



(134) 



dt 



■V^C + AV. 



In the limit e — > 0, with the condition that 



lini -Ko = fim - [p[i] (e) - pp] (e)] 
€—►0 e e^o e 



(135) 



(136) 



exists and is finite, the matrix As r is non-zero only for r — s — 2, s — 1, s, and its eigenvalues an are all negative given 
by 



as = -sAK-i = -s^, s > 1 

where the friction coefficient A — \{y — 0), given by (|62|), is of order e. 
We thus obtain at lowest order in e 



(137) 



1 _ e"^^* 



(138) 



with Vb and A given by Eqs.(|^,^ or <^^. 
At order we find for the case = p[2\ 



{y)t - {y)oo ( 



1 — e 



(139) 
(140) 
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with (V)^ and {V^) ^ given by Eqs.@ |^ with = pp]. 

Similarly for p\^i^ ^ pp] will have another linear combination of exp(— i/ra), exp(— 2i/ra) with 



7\/ + m , , 

Ta = (141) 



We notice that taking realistic numbers one finds Tq of the order 10 ^ — 10 which means that the stationary state 
is reached very rapidly. 



B. Finite systems 



For finite systems the equations for the time evolution of {V)^ and {V"^)^ are given at order by Eqs.( [70|j7l| ) 

where the densities and the coefficients K^^^ are now functions of time. Since we do not consider the full many- 
body problem describing the fluids, we take care of this unknown time dependence by assuming that the fluids satisfy 
the homogeneity condition ([77|). Moreover for , r = —2,-1,1, we take the expressions ( Jt^ ) obtained with the 
Maxwellian distributions. From Eq. ([zs]) we obtain at order e 



di 



2T, 



[1] 




Ira 
M 



TT 



2Tr 



[2] 




8kBT[2] 



(142) 
(143) 



At order e we thus have a system of five O.D.E for the unknown {X)^, {V)^, {V'^')^, T[i] and T[2]. At this order we 
conclude that the system will evolve toward the unique equilibrium state 



r[i](oo) = r[2](oo) = 



N[i] + N[2] 



= L- 



7V[i] + iV[. 


ks ^[i]^[i] 



M 



(144) 

(145) 
(146) 
(147) 



which is identical to the one obtained in thermodynamics for the conducting piston (k 7^ 0). 

Let us then analyze qualitatively how the evolution toward this final equilibrium state takes place. 

Given that at f = 0, we have X = Xp, V = 0, T\^i]{Q) = T^i], T'[2](0) = T[2], and P[i](0) ^P[2](0), then as long as the 
velocity {V)^ remains small Eqs.(142.143) yield the adiabats 



T^2]{t) {L-{X),) 
and the time evolution for the piston is 



T[2] {L ~ X(iY 



^m'^WTZTi^ ^[2] ^[2] 



[L - XqY 

{L-{X)^f 



m 

M 



(148) 
(149) 



(150) 



with A(i) given by Eq.(^, together with Eqs.(|T4|, |l4|). 

The piston evolves thus adiabatically until a time ta, of the order Tq Eq.(141), where the pressures become equal. 



I.e. 



Ar[l]r[i](<,) _ Nl2]T[2]ita) 



L-{X), 



(151) 



which yields 
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w*. = — 

1 + 

T[i]{ta) = T[i] 
T[2]{ta) = 7(2] 



NmT[2] ( L-X„ 
\ 2 



1/3 



' L-Xq 



(152) 

(153) 

(154) 
(155) 



In the next stage, as can be seen from the numerical simulations of Sec. VIII, the pressures will remain approxima- 
tively constant and equal, i.e. 



and 



From Eqs.(156,157) we have 



P[i]{t) = kg ' \' — ^ cte 
P[2]{t) = kg ^^^^ ^ cte 



iV[i]r[i](i) [L-(x)j = iV[2]r[2](<) (X), 



1 d 
1 d 



and using Eqs.( |l42| , |l43| ) we conclude that 

M{v^)^ =fcBy^T[i](<)r[2](i) 



2V?n^ /8fc 
3~M 



In other words the temperatures T\^i^{t) and T[2](t) evolves, but at all time t we have 
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3^ 



(156) 
(157) 

(158) 

(159) 
(160) 

(161) 
(162) 

(163) 
(164) 



where {V)^ and {V'^}^ are the stationary values for the infinite system with temperatures T[,;] — T[i]{t), i — 1,2. 
Introducing the expressions ( |l6l| ) and (162) in ( |l42| ), (143) yields, with = p[2] = p 



and thus 



d 



iT[2] -T[i]) ^ 



2 p /8fc 



3 7V| 



[i]f 



M V TT 



2 p /Sfcs 



I'll] - yTi2] 
T[i] - JT[2] 



dt 



1 \ ^/m /Sfes 



(^[2] - ^[1] 



(165) 
(166) 

(167) 
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Let us note that for = p[2], it follows from Eqs.(|61|,|62|) that 



M ' 



+ A[2] = ^/TO^/■^n[l]T[l 
A 



j.e. K 



M 



■ T\2] 



We thus have 



dt 



{T[2 



Tr 



3 V^[l]fcs iV[2]fcB 



from which follows that the system will reach equilibrium with a relaxation time Tc of the order 

-1 ^ 1 ^ 1 
kB \ ks 

Taking realistic numbers we obtain a relaxation time which is several time the age of the universe. 



(168) 
(169) 
(170) 

(171) 

(172) 



The coefficient 2 



(w[llfc£ 



is clearly the inverse of the specific heat. To understand the origin of the factor 



1/3 we must remember that the piston is moving. Using Eqs.(162.165) yields 

r,[2]^[l] 1 , '^^'[1] , rr \ nllJ-'M 

Pq =2^m''B^+P[i]{V)t = K{T[2]-T[i])=-P^' 

as it should. 



(173) 



Let us also note that the equality of pressure ( 158 ), together with the conservation of energy 

iV[i]T[i](t) + iV[2]r[2](t) + ^ (V)^ = 7V[i]r[i] + iV[2]r[2], 



(174) 



which is exactly satisfied by Eqs.(^,71 142 142) if = p[2], gives the relation between position and temperatures 
(P^,pl7|) 



c (iV[i]T[i] + 7V[2]T[2]) (X), 
L7V[2]r[2](t) ^ (A^[i]T[i] + 7V[2]T[2]) [L - (X),) 



and thus from Eq.(162), or Eq.(|l65|), 



37 Wt = o^\/ Z 



2 /8 /fcB (iV[i]T[i] + iV[2]T[2]) L- {X), 



dt 



3 M V TT 



Finally if we compare ^ (F)^ obtained from Eq.( |l77| ) with Eq.(|70|), we arrive at the conclusion that 



P[2]-P[i]=o(£'(r[2]-T[i])) 



(175) 
(176) 



(177) 



(178) 



which justifies our starting point for the second stage. 



Let us also remark that for some time interval between the adiabatic evolution Eq.(15C), and the approach toward 



equilibrium with p^^] — P[2] (171), there will be an intermediate stage for which Pq 



[l]^[2T 



p, 



[2]-[l] ; 



Q 



is not zero. During 



this intermediate evolution the stochastic motion of the piston will have to be introduced in the thermodynamical 
equations not only with a conductivity coefficient k 7^ 0, but also by means of an internal energy Up and and entropy 
Sp of the piston with Up = Up{Sp). 
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FIG. 2. The average position of the piston as a function of time for the infinite system. The parameters are m — 1, Tjij — 1, 
T[2] = 10, n[i] = 1, n[2] = 1/10 and M = 2, 5, 10, 20, 50, 100 (from top to bottom). The dotted box is enlarged on Fig. 



VIII. NUMERICAL SIMULATIONS 



In order to verify the assumptions on which our previous analysis is based, we made numerical simulations for the 
finite as well as for the infinite system in one dimension taking fc^ = 1. The initial state of the fluids particles of mass 
m is given by Maxwellian distributions of velocities according to Eq. (]l|) . Then the particles and the piston (a particle 
of mass M with initial coordinate {Xq, Vq)) interact only through elastic collisions. 

For the infinite system, we simulate the openness of the system with sources of in-going particles very far from the 
piston position. We compute the average time dependent position {X)t of the piston on 10'^ — 10^ different samples. 
An example of the time evolution of the piston average position is shown on Fig. |^ for the following parameter to = 1, 
T[i] = 1, T[2] — 10, = 1, rt[2] = 1/10, i.e. equal pressure on both sides of the piston, and M = 2, 5, 10, 20 , 50, 100 



(from top to bottom). The initial position and velocity of the piston are Xq = Vq = 0. As expected from Eq.(139) we 
observe that the average position quickly behaves as {X)t = (V)^ {t > 10^). The relaxation time Tq necessary to 
reach this stationary behavior is too short to be represented on Fig. 0, and is presented on Fig. I for M = 5,10,50,100. 
Remark that Tq depends on the ratio m/M . 

Even if the pressures are equal on both sides of the piston, i.e. n[x]T'[i] = '^[2]2^[2]j the stationary velocity (V^)^^ is non 
zero and is presented as a function of M on Fig. ^ for the above parameters, as well as for T\^2] — 20 and n[2] = 1/20. 
The lines on Fig. |^ show the average stationary velocity computed in Sec. VI B up to the fourth order in e = ^ m/M 



A3) 



(179) 



with (see Eq.(102)) 



{V) 



(1) _ 



/27r 1 



(180) 



and 



{V) 



(3) 



27r 1 



V^[l]'^[2] 



(181) 



One sees that the values obtained at the fourth order agree very well with the numerical results. 
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FIG. 3. The average position of thepiston as a function of time for the infinite system. Magnification of the dotted box on 
Fig. |. The parameters are as in Fig. |fwith M = 5, 10, 50, 100. 
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FIG. 4. The stationary average velocity reached by the piston for the infinite system as a function of its mass M. The 
parameters of the simulation are m = 1, T[i] = 1, n[i] = 1 while T[2] = 10, ny2\ = 1/10 (circles) and Tpj = 20, ny2\ = 1/20 
(squares), respectively. The solid, resp. dashed, line shows the theoretical value obtained in the text up to fourth order in e. 
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FIG. 5. The piston's position as a function of time. The parameters of the simulation are m = 1, T[i] — 1, n[i] = 1 while 
Tp] = 10, n[2] = 1/5. X{t) is shown for for M = 2,10,10^,10^,10* (from top to bottom). The dashed line stands for the 
equilibrium solution {X = 3 x 10^) while the dotted lines show the thermodynamical prediction {X — 84260.5) for the adiabatic 
piston {k — 0). 



Concerning the finite system, we present here a qualitative study of the piston properties as a function of the 
parameter e. We consider a finite system of size L = 6 x 10^ with N^i] — N[2] = 10^ and with = 10^ (thus we have 
n[i] = 1, n[2] = 1/5). The initial temperatures on the left and on the right of the piston are Tjij — 1 and T[2] = 10. 
The initial pressure difference is thus p[i] — P[2] = fc-B(^[i]2^[i] — "'[2]r[2]) < 0, pushing the piston toward the left side. 
The time scale of our simulations has to be compared with the typical time it takes for a given particle to cross the 
system which is of the order 0(10^). We plot on Fig. ^ the time behavior of the piston's position for different values 
of e = 2-1/^ 10-"/2 (n = 1, 2, 3, 4). 

On Fig. ^ and Fig. we plotted respectively the evolution of the temperatures Tj^j (t), {i = 1, 2) and of the pressures 
for the same parameters. Notice that the temperatures are estimated through the second moment of the fluid particles 
velocity for all the particles present on the same side of the piston. It is equivalent to the temperatures on the face of the 
piston only if we suppose that the system is homogeneous. The dashed lines show the equilibrium solution predicted 
in SeclvT^ where nfii(oo) = n[2](oo) = (iV[i] +7V[2])/L, = L7V[i]/(iV[i] +iV[2]), T[i](oo) = r[2](oo) = (T^i] +r[2])/2, 
P[i](oo) = p[2](oo) = kB{N[i] +iV[2])(T[i] +T[2])/(2L). On the other hand, the dotted lines show the equilibrium values 
obtained from a numerical integration of the phenomenological thermodynamical equations ( |39| , ^,p3|) with k = 
for the primitive model of an adiabatic piston and where we used A[i](y) ~ \i]{V = 0) = n^q {t)^y8kBT^i]m/n. In 
this case the system evolves toward a metastable state with p[^{oo) = p[2](c>o) but T[^{oo) ^ T^2]ioo). One should 
note that on the scale of Figs. |5[ ^ and 0, one can hardly distinguish the results obtained by numerical integr ation 



84260.5, Tfi] ~ 1.545, T[2] ~ 9.455 and = pp] 



1.83) from those obtained by heuristic arguments Eqs.(152 
jl55| ), X ~ 82200, T[l] ~ 1.48, T^] ~ 9.32 and = pp] ~ 1.80. 
Let us first comment on the Fig. ^ where one sees that the pressures quickly become constant and equal to their 
final values. The time scale ta needed to equalize the pressures is roughly independent on the mass ratio m/M . This 
behavior has been used in Sec. VII B| to give a qualitative discussion of the evolution. Consider now the Fig. |^ where 
one sees that the temperatures reach quickly an intermediate stage characterized by the equality of the pressures (time 
scale to) then evolve slowly toward their equilibrium values on a time scale t^. The time scale To depends strongly on 
the mass ratio = m/M and for e <§C 1 we have Te 3> ta- One remarks that, for times covered in our simulations and 
for ~ 10^'' the system behaves as if the piston was "truly" adiabatic (k = 0). 

We plotted on Fig. || and Fig. ^ the velocity distribution for the particles of both the fluids. It is averaged on 

21 = 10, P\l] = P[2\ = 1, 

We present 



all particles located on the same side of the piston. The initial conditions are T[i] = \, T\ 



= Vb = 0, m = 1 and M — 10. The dotted lines show the Maxwellian initial distribution 
these distributions for three different times t = ti,t2,t3. 



10, P[i] 
6T[.j iv). 



24 



0.0 I — 
Oe+00 



le+06 
t 



2e+06 



FIG. 6. The temperature of the fluids on both sides of the piston as a function of time. The parameters of the simulation 
are m = 1, T[i] = 1, ny^ = 1 while Tp] = 10, n[2] = 1/5. 7(1] (i), resp. Ty2]{t), is shown for for M = 2, 10, 10^ 10'\ lO" (from 
middle to bottom, resp. to top). The dashed line stands for the equilibrium solution (T[i](oo) = r[2j(oo) = (T[i] +T[2])/2 = 5.5) 
while the dotted lines show the thermodynamical prediction (T[i](oo) ~ 1.545, T[2](oo) ^ 9.455) for k = 0. 




FIG. 7. The pressures of the fluids on both sides of the piston as a function of time. The parameters of the simulation 

are m = 1, T[i] = 1, n[i] = 1, = 1 while Tpj = 10, fi[2] = 1/5, ppj = 2. P[i](t), resp. P[2](t), is shown for for 
M = 10, 10^, 10^, lO**. The dashed line stands for the equilibrium solution which is equal to the thermodynamical prediction 
(p[i](oo) = p[2](oo) = iV[i](T[i] + Tp])/L = 11/6). 
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FIG. 8. The particles velocity distribution <j)[i]{v,t) of the fluid [1] at different time for the finite system. The initial 
conditions arc Tfij = 1, T[2] = 10, P[i] = p[2] ~ 1, = Vb = 0, m = 1 and M = 10. The dotted line show the initial 
distribution (pTf^i] while the dashed line show the expected equilibrium distribution ipTi^ (oo) {v) with Tfij (oo) = (r[i] + T[2] ) /2. 
The squares show the velocity distribution for ti = 2 x 10^ where the temperature of the fluid is T[i]{ti) ~ 1.6. The solid line 
is for t2 = 2.5 X 10® with T[i](t2) ^ 3.7 while the circles are for ta = 7.5 x 10^ where T[i](t3) ^ T'[i](oo)- 




FIG. 9. The particles velocity distribution (j>[2]{v,t) of the fluid [2] at different time for the finite system. The initial 
conditions are r[i] = 1, T[2j = 10, = p[2] = 1, Xo = Vo = 0, m = 1 and M = 10. The dotted line show the initial 
distribution (^Tp] (v) while the dashed line show the expected equilibrium distribution (f>^2]{oo){v) with T[2](oo) = (T[i] +T[2])/2. 
The squares show the velocity distribution for ii = 2 x 10^ where the temperature of the fiuid is T[2]{ti) ~ 9.4. The solid line 
is for t2 = 2.5 X 10'^ with T[2](t2) ^ 7.3 while the circles are for ts = 7.5 x 10^ where TpjCis) ^ r[2](oo). 
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FIG. 10. The time evolution of the temperatures T'[i](t) (bottom) and T[2]{t) (up) for the same initial conditions as for Figs. 
^and|^. The successive vertical dotted lines show the times ti,t2 and at which the velocity distributions plotted on Figs. ^ 
and ^ were determined. 

For the time ti = 2 x 10^ the system has not reached the relaxation time ta and the pressures are not yet constant. 
At the time t2 {ta < t2 < Tq), the system has constant pressures but has not reached the equilibrium state. At 
time is > To the system has reached its final equilibrium state. We present on Fig. |l^ the time dependence of the 
fluid temperatures. Vertical lines stand for times ii,i2,^3- We see that, for times such as t < Tc, these distributions 
are not Maxwellian but develops peaks or wells for small velocities. Moreover these intermediate time distribution 
are clearly not symmetric. This facts are the reason why we did not suppose Maxwellian distributions for the fluid 
particles in our previous analytical study. After a sufficiently long time of the order of the equilibrium relaxation 
time of the system Tg, the distributions are very well fitted with the Maxwellian distribution and wc conjecture that 
(j>[i]{v,oo) = 0[2](w,oo) = 0T[i,(oo)(^^) with T[i](oo) = (T[i] +T[2])/2 = 11/2. The dashed lines on Fig. |and Fig. | 
show the conjectured equilibrium distribution. This last conjecture is valid only for m < M as it has been shown that 
the equilibrium distribution in the case of m = M is a simple superposition of the two initial Maxwellian j^] . 

IX. CONCLUSIONS 

Several conclusions can be drawn from this investigation on the stochastic motion of a piston which is adiabatic 
when rigidly fixed. 

For infinite systems, it has been shown that very rapidly a stationary state is reached, where the average velocity of 
the piston is a function of the pressures and temperatures of the fluids on both sides. In particular if the pressures are 
equal but the temperatures different the final average velocity is non-zero and the piston has a macroscopic motion 
toward the high temperature region. This result is related to the asymmetry in the fluctuations of the piston, due 
to the asymmetry of the temperatures on both sides, which in turn induces both a macroscopic force and a heat 
current from one side to the other. In other words in the stationary state the stochastic piston is a conductor. To 
obtain a stationary state with zero average velocity it is necessary to compensate the stochastic force by a macroscopic 
difference in the pressures of the fluids. In this very special state no work is done by the piston on the fluids, but 
heat is transferred from one side to the other. Moreover this state is not a state of mechanical equilibrium since all 
the odd moments of the piston velocity, except the first one, are non zero. 

For finite systems it has been shown that the system will always evolve toward the equilibrium state predicted for 
a conducting piston where pressures, temperatures and densities are equal on both sides of the piston , but the time 
needed to reach this final equilibrium state will depend strongly on the mass ratio m/M and can reach several time 
the age of the universe for reasonable numbers. This evolution takes place in two stages. In the first stage the system 



27 



reaches rather rapidly (although not as fast as the time needed to reach the stationary state of the infinite system) a 
state of "metastable equilibrium" where the pressures are equal but the temperatures different. This initial evolution 
corresponds to an adiabatic evolution (no heat transfer) and the results obtained from the microscopical theory, from 
the thermodynamic equations, and from the numerical simulations are in good agreement. In the second stage the 
piston and the temperatures of the fluids evolve in such a way that the pressures remain approximatively constant 
and equal. If m <C M it appears that at all time the average velocity of the piston coincide with the velocity in the 
stationary state of the infinite system with values of temperatures and pressures given by those of the finite system 
at that time. This observation can be understood from the fact that the stationary state of the infinite system is 
reached very rapidly. If m, <C M the time needed to reach the final state is so large that on the numerical simulations 
the "metastable equilibrium" state appear to be stable and the piston behaves as an adiabatic piston on the time 
scale considered. On the other hand for m ~ M the final equilibrium state is reached rather rapidly. Moreover it 
was observed that during the time evolution the velocity distribution of the fluid particles is not Maxwellian and 
that this distribution evolves slowly toward the equilibrium Maxwellian distribution characterized by the equilibrium 
temperature. 

Finally the conclusions obtained from thermodynamics, from kinetic theory, and from numerical simulations, are 
all in good agreements. There appears to be no violation of the second law and no paradox involved. 
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